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Abstract 

The inclusive production of D* ± mesons in photon-photon collisions has been measured using 
the OPAL detector at LEP at e + e~ centre-of-mass energies \/s ee of 183 and 189 GeV. The D* + 
mesons are reconstructed in their decay to D°7r + with the D° observed in the two decay modes 
K~7r + and K~7r + 7r~7r + . After background subtraction, 100.4 ± 12.6 (stat) D* ± mesons have 
been selected in events without observed scattered beam electron ("anti-tagged") and 29.8 ± 
5.9 (stat) D* 1 * 1 mesons in events where one beam electron is scattered into the detector ("single- 
tagged"). Direct and single- resolved events are studied separately. Differential cross-sections 
dcr/dpT and d<r/d|?7 D | as functions of the D* ± transverse momentum pijT an d pseudorapidity 
rp* are presented in the kinematic region 2 GeV < < 12 GeV and |?7 D | < 1.5. They 
are compared to next-to-leading order (NLO) perturbative QCD calculations. The total cross- 
section for the process e + e~— > e + e~cc where the charm quarks are produced in the collision of 
two quasi-real photons is measured to be a(e + e~ — > e + e~cc) = 842 ± 97 (stat) ± 75 (sys) ± 
196 (extrapolation) pb. A first measurement of the charm structure function F^ c of the photon 
is performed in the kinematic range 0.0014 < x < 0.87 and 5 GeV 2 < Q 2 < 100 GeV 2 , and the 
result is compared to a NLO perturbative QCD calculation. 



(To be submitted to Eur. Phys. J. C) 



The OPAL Collaboration 



G. Abbiendi 2 , K. Acker staff 8 , P.F. Akesson 3 , G.Alexander 23 , J.Allison 16 , K.J.Anderson 9 , 
S.Arcelli 17 , S.Asai 24 , S.F.Ashby 1 , D. Axen 29 , G. Azuelos 18 ' a , I.Bailey 28 , A.H.Ball 8 , 
E. Barberio 8 , T. Barillari 2 , R.J. Barlow 16 , J.R. Batley 5 , S.Baumann 3 , T.Behnke 27 , K.W. Bell 20 , 
G.Bella 23 , A. Bellerive 9 , S. Bentvelsen 8 , S. Bethke 14 ^, S.Betts 15 , O. Biebel 14 ^, A.Biguzzi 5 , 
I. J. Bloodworth 1 , P. Bock 11 , J. B6hme 14 ' h , O.Boeriu 10 , D. Bonacorsi 2 , M. Boutemeur 33 , 
S.Braibant 8 , P. Bright-Thomas 1 , L. Brigliadori 2 , R.M. Brown 20 , H. J. Burckhart 8 , 
PCapiluppi 2 , R.K.Carnegie 6 , A.A. Carter 13 , J.R. Carter 5 , C.Y.Chang 17 , D.G. Charlton 1,6 , 

D. Chrisman 4 , C. Ciocca 2 , P.E.L. Clarke 15 , E.Clay 15 , I.Cohen 23 , J.E. Conboy 15 , O.C.Cooke 8 , 
J.Couchman 15 , C. Couyoumtzelis 13 , R.L.Coxe 9 , M.Cuffiani 2 , S.Dado 22 , G.M. Dallavalle 2 , 
S.Dallison 16 , R.Davis 30 , A. de Roeck 8 , PDervan 15 , K. Desch 27 , B. Dienes 32,h , M.S. Dixit 7 , 

M.Donkers 6 , J.Dubbert 33 , E. Duchovni 26 , G.Duckeck 33 , LP Duerdoth 16 , PG. Estabrooks 6 , 

E. Etzion 23 , F.Fabbri 2 , A. Fanfani 2 , M.Fanti 2 , A.A.Faust 30 , L.Feld 10 , P.Ferrari 12 , 

F. Fiedler 27 , M. Fierro 2 , I. Fleck 10 , A.Frey 8 , A.Fiirtjes 8 , D.I. Futyan 16 , P. Gagnon 12 , 

J.W.Gary 4 , G. Gaycken 27 , C. Geich-Gimbel 3 , G. Giacomelli 2 , P. Giacomelli 2 , 
D.M. Gingrich 30 ' , D. Glenzinski 9 , J. Goldberg 22 , W. Gorn 4 , C. Grandi 2 , K. Graham 28 , 

E. Gross 26 , J. Grunhaus 23 , M. Gruwe 27 , C.Hajdu 31 G.G.Hanson 12 , M.Hansroul 8 , M.Hapke 13 , 
K. Harder 27 , A.Harel 22 , C.K.Hargrove 7 , M. Harin-Dirac 4 , M. Hauschild 8 , C.M.Hawkes 1 , 
R. Hawkings 27 , R.J.Hemingway 6 , G.Herten 10 , R.D. Heuer 27 , M.D. Hildreth 8 , J.C.Hill 5 , 

P.R. Hobson 25 , A.Hocker 9 , K.Hoffman 8 , R.J. Homer 1 , A.K.Honma 8 , D. Horvath 31,c , 
K.R. Hossain 30 , R.Howard 29 , P. Hiintemeyer 27 , P. Igo-Kemenes 11 , D.C.Imrie 25 , K.Ishii 24 , 

F. R. Jacob 20 , A. Jawahery 17 , H. Jeremie 18 , M. Jimack 1 , C.R. Jones 5 , P. Jovanovic 1 , T.R. Junk 6 , 
N.Kanaya 24 , J. Kanzaki 24 , G. Karapetian 18 , D.Karlen 6 , V. Kartvelishvili 16 , K.Kawagoe 24 , 
T.Kawamoto 24 , P.I.Kayal 30 , R.K. Keeler 28 , R.G. Kellogg 17 , B.W.Kennedy 20 , D.H.Kim 19 , 

A.Klier 26 , T. Kobayashi 24 , M.Kobel 3 , T.P.Kokott 3 , M.Kolrep 10 , S. Komamiya 24 , 
R.V. Kowalewski 28 , T.Kress 4 , P. Krieger 6 , J. von Krogh 11 , T.Kuhl 3 , M. Kupper 26 , P. Kyberd 13 , 

G. D. Lafferty 16 , H. Landsman 22 , D. Lanske 14 , J. Lauber 15 , I. Lawson 28 , J.G.Layter 4 , 
D.Lellouch 26 , J.Letts 12 , L. Levinson 26 , R. Liebisch 11 , J.Lillich 10 , B.List 8 , C. Littlewood 5 , 

A.W.Lloyd 1 , S.L.Lloyd 13 , F.K. Loebinger 16 , G.D.Long 28 , M. J. Losty 7 , J.Lu 29 , J. Ludwig 10 , 
A. Macchiolo 18 , A. Macpherson 30 , W. Mader 3 , M.Mannelli 8 , S. Marcellini 2 , T.E. Marchant 16 , 
A.J.Martin 13 , J.P.Martin 18 , G.Martinez 17 , T. Mashimo 24 , P. Mattig 26 , W.J. McDonald 30 , 
J.McKenna 29 , E.A. Mckigney 15 , T.J. McMahon 1 , R.A. McPherson 28 , F.Meijers 8 , 
P.Mendez-Lorenzo 33 , F.S. Merritt 9 , H. Mes 7 , I.Meyer 5 , A. Michelini 2 , S.Mihara 24 , 

G. Mikenberg 26 , D.J.Miller 15 , W.Mohr 10 , A. Montanari 2 , T.Mori 24 , K.Nagai 8 , I. Nakamura 24 , 

H. A. Neal 12 ^, R. Nisius 8 , S.W. O'Neale 1 , F.G. Oakham 7 , F. Odorici 2 , H.O. Ogren 12 , 

A. Okpara 11 , M.J. Oreglia 9 , S. Orito 24 , G.Pasztor 31 , J.R. Pater 16 , G.N.Patrick 20 , J. Patt 10 , 
R.Perez-Ochoa 8 , S. Petzold 27 , P. Pfeifenschneider 14 , J.E. Pilcher 9 , J. Pinfold 30 , D.E. Plane 8 , 

B.Poli 2 , J.Polok 8 , M. Przybycieh 8,d , A. Quadt 8 , C.Rembser 8 , H.Rick 8 , S.A.Robins 22 , 
N.Rodning 30 , J.M. Roney 28 , S.Rosati 3 , K. Roscoe 16 , A.M.Rossi 2 , Y. Rozen 22 , K.Runge 10 , 
O.Runolfsson 8 , D.R.Rust 12 , K.Sachs 10 , T.Saeki 24 , O.Sahr 33 , W.M.Sang 25 , 
E.K.G.Sarkisyan 23 , C.Sbarra 28 , A.D. Schaile 33 , O.Schaile 33 , P. Scharff-Hansen 8 , J.Schieck 11 , 
S. Schmitt 11 , A. Schoning 8 , M. Schroder 8 , M. Schumacher 3 , C. Schwick 8 , W.G. Scott 20 , 
R. Seuster 14 ' fe , T.G. Shears 8 , B.C. Shen 4 , C.H. Shepherd-Themistocleous 5 , P.Sherwood 15 , 
G.P Siroli 2 , A. Skuja 17 , A.M. Smith 8 , G.A. Snow 17 , R. Sobie 28 , S. Soldner-Rembold 10 ' 6 , 
S. Spagnolo 20 , M. Sproston 20 , A. Stahl 3 , K. Stephens 16 , K. Stoll 10 , D. Strom 19 , R. Strohmer 33 , 

B. Surrow 8 , S.D.Talbot 1 , P. Taras 18 , S.Tarem 22 , R. Teuscher 9 , M. Thiergen 10 , J.Thomas 15 , 



I 



M.A.Thomson 8 , E. Torrence 8 , S. Towers 6 , T. Trefzger 33 , I. Trigger 18 , Z. Trocsanyi 32 ' 9 , 
E. Tsur 23 , M.F. Turner- Watson 1 , I.Ueda 24 , R. Van Kooten 12 , P. Vannerem 10 , M. Verzocchi 8 , 
H.Voss 3 , F.Wackerle 10 , D.Waller 6 , C.P.Ward 5 , D.R.Ward 5 , P.M. Watkins 1 , A.T.Watson 1 , 
N.K.Watson 1 , P.S.Wells 8 , T.Wengler 8 , N.Wermes 3 , D. Wetterling 11 J.S.White 6 , 

G.W.Wilson 16 , J.A.Wilson 1 , T.R.Wyatt 16 , S. Yamashita 24 , V. Zacek 18 , D.Zer-Zion 8 



1 School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 
2 Dipartimento di Fisica dell' Universita di Bologna and INFN, 1-40126 Bologna, Italy 
3 Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 
4 Department of Physics, University of California, Riverside CA 92521, USA 
5 Cavendish Laboratory, Cambridge CB3 OHE, UK 

6 Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ottawa, 
Ontario K1S 5B6, Canada 

7 Centre for Research in Particle Physics, Carleton University, Ottawa, Ontario K1S 5B6, 
Canada 

8 CERN, European Organisation for Particle Physics, CH-1211 Geneva 23, Switzerland 
9 Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 60637, 
USA 

10 Fakultat fur Physik, Albert-Ludwigs-Universitat, D- 79104 Freiburg, Germany 

11 Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 

12 Indiana University, Department of Physics, Swain Hall West 117, Bloomington IN 47405, 

USA 

13 Queen Mary and Westfield College, University of London, London El 4NS, UK 
14 Technische Hochschule Aachen, III Physikalisches Institut, Sommerfeldstrasse 26-28, D-52056 
Aachen, Germany 

15 University College London, London WC1E 6BT, UK 

16 Department of Physics, Schuster Laboratory, The University, Manchester M13 9PL, UK 
17 Department of Physics, University of Maryland, College Park, MD 20742, USA 
18 Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, 
Canada 

19 University of Oregon, Department of Physics, Eugene OR 97403, USA 
20 CLRC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK 
22 Department of Physics, Technion- Israel Institute of Technology, Haifa 32000, Israel 
23 Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
24 International Centre for Elementary Particle Physics and Department of Physics, University 
of Tokyo, Tokyo 113-0033, and Kobe University, Kobe 657-8501, Japan 

25 Institute of Physical and Environmental Sciences, Brunei University, Uxbridge, Middlesex 
UB8 3PH, UK 

26 Particle Physics Department, Weizmann Institute of Science, Rehovot 76100, Israel 
27 Universitat Hamburg/DESY, II Institut fur Experimentalphysik, Notkestrasse 85, D-22607 
Hamburg, Germany 

28 University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, Canada 
29 University of British Columbia, Department of Physics, Vancouver BC V6T 1Z1, Canada 
30 University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 
31 Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, Hungary 
32 Institute of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 



2 



33 Ludwigs-Maximilians-Universitat Munchen, Sektion Physik, Am Coulombwall 1, D-85748 
Garching, Germany 

a and at TRIUMF, Vancouver, Canada V6T 2A3 

6 and Royal Society University Research Fellow 

c and Institute of Nuclear Research, Debrecen, Hungary 

d and University of Mining and Metallurgy, Cracow 

e and Heisenberg Fellow 

' now at Yale University, Dept of Physics, New Haven, USA 

9 and Department of Experimental Physics, Lajos Kossuth University, Debrecen, Hungary 
h and MPI Munchen 

1 now at MPI fur Physik, 80805 Munchen. 



1 Introduction 

Charged D* mesons^ provide a clean tag to study open charm production in photon-photon 
collisions. The inclusive cross-section for the production of D* mesons can be calculated in 
perturbative QCD (pQCD). Since the process is characterised by two distinct scales, the mass 
m c and the transverse momentum p^ of the charm quarks, two different approaches exist for the 
next-to-leading order (NLO) pQCD calculations. In the "massless" scheme, charm is treated 
as an active flavour in the photon. This approach is expected to be valid for p^ ^> m c . In the 
"massive" scheme, matrix elements for massive charm quarks are used and no charm content 
is assigned to the parton distributions of the photon. One expects this ansatz to be valid at 
Pt ~ rn c [|TJ. 

A third scale is the four-momentum squared, Qf =l2 , of the interacting virtual photons. In 
this paper, two kinematic cases are studied, depending on Qf. For the largest part of the 
cross-section, both exchanged photons are quasi- real (Q\,Q\ ~ 0) and the beam electrons are 
scattered at very small angles. These events are selected by rejecting events with a scattered 
electron in the detector. Events of this type are called anti-tagged. If a photon is highly virtual 
{Ql ^> 0), the corresponding beam electron is usually scattered into the acceptance of the 
detector. Events with one detected scattered electron are called single-tagged. 

In direct events, the two photons couple directly to the cc pair. In resolved events, one pho- 
ton ( "single- resolved" ) or both photons ( "double- resolved" ) fluctuate into a hadronic state and 
a gluon or a quark of the hadronic fluctuation of the photon takes part in the hard interaction. 
For anti-tagged events at e + e~ centre-of-mass energies -y/s ee around 189 GeV, the production of 
D* mesons in photon-photon collisions in leading order (LO) QCD proceeds mainly via direct 
(77 ~~ > cc) and single-resolved (g7 — > cc) photon-photon processes ||[3|], whereas the contri- 
bution from double-resolved events (gg — > cc) is expected to be small. The measurement of 
the open charm cross-section is therefore expected to be sensitive to the gluon content of the 
photon through the photon-gluon fusion process g7 — > cc. The production of D* mesons from 
open bottom production in photon-photon events is expected to be suppressed by more than 

throughout this paper D* refers to D* + as well as to D*~. Charge conjugated modes are always implied. 
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two orders of magnitude 0. Bottom production is suppressed due to the smaller electric charge 
and the larger mass of the b quarks. 



The single-tagged process can be regarded as deep inelastic electron-photon scattering. In 
this configuration, the electron radiating the highly virtual photon, 7*, probes the structure of 
the quasi-real photon, 7, radiated from the second beam electron, and allows the determination 
of the photon structure function F^(x, Q 2 ). Here x is the Bjorken scaling variable and Q 2 the 
virtuality of the highly virtual photon. Due to the large mass of the charm quark, the charm 
structure function of the photon can be calculated in pQCD to NLO ||. In QCD, F% c 
receives contributions from the point-like and the hadron-like structure of the quasi-real photon. 
The two contributions are expected to be well separated in x, with the point-like contribution 
dominating at high x and the hadron-like component sizeable only for x < 0.1. 

The production of D* mesons in photon-photon events has been measured previously by 
JADE §, TASSO § , TPC/27 0, TOPAZ f§, AMY @, ALEPH and L3 H at y/s ee 



ranging from 29 GeV up to 189 GeV. The charm structure function F^ c of the photon has never 
been measured before. The analysis presented here uses the data taken with the OPAL detector 
in 1997 at \/s ee = 183 GeV and in 1998 at i/i ee = 189 GeV with integrated luminosities £ of 
55 and 165 pb -1 , respectively. 



2 The OPAL detector 



A detailed description of the OPAL detector can be found in Ref. [12], and therefore only a 



brief account of the main features relevant to the present analysis will be given here. 

The central tracking system is located inside a solenoidal magnet which provides a uniform 
axial magnetic field of 0.435 T along the beam axisQ. The magnet is surrounded by a lead-glass 
electromagnetic calorimeter (ECAL) and a hadronic sampling calorimeter (HCAL). Outside 
the HCAL, the detector is surrounded by muon chambers. There are similar layers of detectors 
in the endcaps. The region around the beam pipe on both sides of the detector is covered by 
the forward calorimeters and the silicon-tungsten luminometers. 

Starting with the innermost components, the tracking system consists of a high precision 
silicon microvertex detector (SI), a precision vertex drift chamber (CV), a large volume jet 
chamber (CJ) with 159 layers of axial anode wires and a set of z chambers measuring the track 
coordinates along the beam direction. The transverse momenta of tracks are measured with a 



precision of <t Pt /pt = y0.02 2 + (0.0015 ■ Pt) 2 G°t in GeV)0in the central region | cos^l < 0.73. 
In this paper, transverse is always defined with respect to the z direction of the detector. The 
jet chamber also provides energy loss measurements which are used for particle identification. 

The ECAL completely covers the azimuthal range for polar angles satisfying | cos#| < 0.98. 
The barrel section, which covers the polar angle range | cos^l < 0.82, consists of a cylindrical 



2 In the OPAL coordinate system the x axis points towards the centre of the LEP ring, the y axis points 
upwards and the z axis points in the direction of the electron beam. The polar angle 6 is defined with respect 
to the z axis. 

3 Throughout this paper we use the convention c = 1. 
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array of 9440 lead-glass blocks with a depth of 24.6 radiation lengths. The endcap sections 
consist of 1132 lead-glass blocks with a depth of more than 22 radiation lengths, covering the 
polar angle between 0.81 < | cos^l < 0.98. 

The forward calorimeters (FD) at each end of the OPAL detector consist of cylindrical 
lead-scintillator calorimeters with a depth of 24 radiation lengths divided azimuthally into 16 
segments. The electromagnetic energy resolution is about 18%/\^E~ (E in GeV). The acceptance 
of the forward calorimeters covers the angular range from 47 to 140 mrad from the beam 
direction. Three planes of proportional tube chambers at 4 radiation lengths depth in the 
calorimeter measure the directions of electron showers with a precision of approximately 1 mrad. 

The silicon tungsten detectors (SW) |T3| at each end of the OPAL detector cover an angular 
region between 33 and 59 mrad in front of the forward calorimeters. Each calorimeter consists of 
19 layers of silicon detectors and 18 layers of tungsten, corresponding to a total of 22 radiation 
lengths. Each silicon layer consists of 16 wedge-shaped silicon detectors. The electromagnetic 
energy resolution is about 25%/\^E (E in GeV). The radial position of electron showers in the 
SW calorimeter can be determined with a typical resolution of 0.06 mrad in the polar angle 9. 



3 Process kinematics 



The kinematic properties of the two interacting photons are described by their negative squared 
four- momentum transfers, Q 2 =12 i which are related to the scattering angles 6[ relative to the 
beam direction of the corresponding electrons by 

Ql = -{h - k[f « 2^(1 - cos 01), (1) 

neglecting the mass m e of the electron. The quantities ki and k[ are the four-momenta of 
the beam and scattered electrons, and E^ and E[ are their respective energies. The flux of 
transversely polarized quasi-real photons with an energy fraction z of the beam energy and a 
negative squared four- momentum, denoted with Q 2 , may be obtained by the Equivalent Photon 
Approximation (EPA) 0: 

d 2 iV 7 a_ ( l + (l-z) 2 _ 2m|z\ 

dzdQ 2 ~ 2n{ zQ 2 Qf )' U 

where a is the fine structure constant. The minimum kinematically allowed squared four- 
momentum transfer, Q min , is determined by the electron mass, 

2 Tfl 2 Z 2 

Qmin = -y _ • (3) 

The maximum squared four- momentum transfer, Q max , is given by the experimental anti- 
tagging condition according to Eq. [1]. 

If one of the photons is highly virtual (Q\ = Q 2 = —q 2 > Ql = P 2 = —p 2 « 0), the 
cross-section of the process e + e~ — > e + e~cc can be written as a product of the deep inelastic 
electron-photon scattering cross-section and the flux of quasi-real photons, 

d (Tg+g-^g+g-^ d (jg^^ggg d A^y 



dxdg 2 d^dP 2 dxdQ 2 d^dP 2 ' 
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The differential deep inelastic electron-photon scattering cross-section may be parametrised in 
terms of structure functions as JL5| 



dxdQ 2 xQ 4 



(1 + (1 - yf)FUx, Q 2 ) - y 2 FUx, Q 2 ) , (5) 



where x = Q 2 /2pq and the inelasticity y = pq/pki are the usual dimensionless deep inelastic 
scattering variables. In the kinematic regime studied in this paper, ?/ 2 <C 1. The contribution 
proportional to F^ c (x,Q 2 ) in Eq. ^ is therefore neglected. The variable x is experimentally 
accessible via the relation 

= Q2 ~ Q2 

X Q 2 + W 2 + P 2 ~ Q 2 + W 2 ' { ) 

where W 2 is the invariant mass squared of the photon-photon system. The event selection 
ensures that the virtuality P 2 is usually very small compared to Q 2 , so P 2 is neglected for 
the determination of x from Eq. |6|. By measuring the deep inelastic electron-photon scattering 
cross-section for e + e~ — > e + e" CC cLS £L function of x and Q 2 , the charm structure function 
F2 C (x, Q 2 ) of the photon can be determined. 



4 Monte Carlo simulation 



For real photons the PYTHIA 6.121 jl6] photon-photon Monte Carlo program, based on LO 
pQCD calculations, is used to simulate the process e + e~— > e + e~77 — > e + e~cc — > e + e~D*X 
(X is a hadronic system). Two distinct samples, one for the direct process, 77 — ► cc, and 
one for the single-resolved process, g7 — > cc, were generated using matrix elements for massive 
charm quarks. The different e + e~ centre-of-mass energies, \fs ee = 183 and 189 GeV, are taken 
into account by generating events at both energies according to the ratio of the corresponding 
integrated luminosities. In case of the single-resolved process, the SaS-lD parametrisation [|Tjj 



is used for the parton distributions of the photon. The fragmentation of the charm quarks is 
modelled using the Peterson fragmentation function |T8| with the PYTHIA default parameter 



e c = 0.031, and the charm mass is set to m c = 1.6 GeV. The resulting average scaled energy 
of the generated D* mesons is (xd*) = 2 ■ (Ed* /W) = 0.84, where Ed* is the D* energy in 
the photon-photon centre-of-mass system. Final state QCD radiation off the primary charm 
quarks is taken into account using the leading-log approximation. A sample of double-resolved 
quasi-real photon-photon events (gg — > cc) was also generated with the PYTHIA 6.121 Monte 
Carlo generator. 



The LO Monte Carlo generators HERWIG 5.9 and Vermaseren |20| are used to model 
the D* production in deep inelastic electron-photon scattering, e + e~— > e + e~7*7 — > e + e~cc — > 
e + e~D*X. For both Monte Carlo generators, the charm quark mass is chosen to be m c = 
1.5 GeV. In HERWIG, the cross-section is evaluated for massless charm quarks. The charm 
production is modelled using matrix elements for massless charm quarks, together with the GRV 
parametrisation [^1] for the parton distributions of the photon, again for massless charm quarks. 
The effect of the charm quark mass is accounted for rather crudely by not simulating events 
with W < 2m c , giving an unphysically sharp step in the cross-section at this threshold. Due to 
the massless approach used in HERWIG and the crude treatment at threshold, the predicted 
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charm production cross-section is likely to be too large. Nevertheless, the final state kinematics 
are treated correctly, so HERWIG can be used for the determination of selection efficiencies. 
The fragmentation of quarks into hadrons is modelled via the cluster fragmentation model 
yielding an average scaled energy of the generated D* mesons of (xd*) = 0.64. The Vermaseren 
generator is based on the Quark Parton Model (QPM) and consequently does not take into 
account the hadron-like component of the photon structure. It models the complete dependence 
of the cross-section on the different photon helicities. The fragmentation into hadrons is handled 
via JETSET 7.4 |Uj, where the same fragmentation model is used as for the PYTHIA Monte 
Carlo described above. The average scaled energy of the generated D* mesons is (xd*) — 0.82, 
close to the PYTHIA value. 

The e + e~ annihilation background with D* mesons in the final state has been simulated with 
the PYTHIA 5.7 Monte Carlo model. The Monte Carlo generator GRC4F |22| was used to 



simulate four-fermion processes that are background to the photon-photon sample. All Monte 
Carlo samples were generated with full simulation of the OPAL detector [^3j. They are analysed 
using the same reconstruction algorithms as applied to the data. 



5 Event selection and D* reconstruction 



In this and in the following three sections, only anti-tagged photon-photon scattering events 
are studied; the analysis of tagged electron-photon scattering events is described in section 9. 
Because the difference of the e + e~ centre-of-mass energies of \fs ee = 183 and 189 GeV is small 
the data samples recorded at both energies are combined. Anti-tagged photon-photon events 
are selected using the following set of cuts: 

• At least three tracks must have been found in the tracking chambers (SI, CV and CJ). 
A track is required to have a minimum transverse momentum of 120 MeV, more than 20 
hits in the central jet chamber used to calculate the specific energy loss dE'/dx, and the 
innermost hit of the track must be inside a radius of 60 cm with respect to the z axis. 
The distance of closest approach to the origin must be less than 20 cm in the z direction 
and less than 1 cm in the r<j) plane. 

• To reduce background from e + e~ annihilation events with D* mesons in the final state, the 
sum of all energy deposits in the ECAL is required to be less than 40 GeV. Calorimeter 
clusters have to pass an energy threshold of 100 MeV for the barrel section and 250 MeV 
for the endcap sections. 

• To reduce the e + e~ annihilation background further, the visible invariant mass of the 
event, W v i s , should be less than 60 GeV. W vis is calculated using the energies and positions 
of clusters measured in the ECAL, the HCAL, the FD and the SW calorimeters and using 
the momenta of tracks. A matching algorithm is applied to avoid double-counting of 
particle momenta in the calorimeters and in the tracking chambers. 

• Anti-tagged events are selected by vetoing all events containing an energy deposit of more 
than 50 GeV in the SW or FD in either hemisphere of the detector. This corresponds to 
a maximum allowed scattering angle of the beam electrons of 9' = 33 mrad for electrons 
with E' > 50 GeV. 
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The method of reconstructing D* mesons is similar to that used in former OPAL analy- 
ses [25]. It exploits the small mass difference between the D* and the D° mesons which causes 
the kinetic energy of the slow pion in the decay D* + — > D°7r + to be only 6 MeV in the D* rest 
frame. Thus, the combinatorial background is small due to the limited phase space. The D° 
mesons are identified via their decay D° -> K"7r + and D° -> K" ~7r + 7r 7r + , which form, together 
with the slow pion, the "3-prong" and "5-prong" decay modes of the D*, respectively. 

In the 3-prong decay mode, all combinations of two oppositely charged tracks in an event 
are used to form D° candidates. The dE/dx probability j>^^ dx f or the kaon hypothesis should 
exceed 10% for at least one of the two tracks. The invariant mass M^o nd of this combination 
is calculated, assigning the kaon mass to the kaon candidate and the pion mass to the other 
track. If for both tracks , p^/ dx j s greater than 10%, both possible Kn combinations are used. 

If M^J Ld lies within a window around the nominal D° mass, 

1790 MeV < M^ a nd < 1940 MeV, (7) 

the combination is retained as a D° candidate. All remaining tracks of opposite charge to 
the kaon candidate are then examined and the invariant mass Mf^ nd of the D* candidates is 
calculated assigning the pion mass to the third track. 

Random combinations of low- momentum tracks are the largest source of background passing 
the above cuts. To reduce this background, we exploit the fact that the D° is a pseudo-scalar 
particle which decays isotropically in its rest frame. This leads to a flat distribution of cos#*, 
where 9* denotes the decay angle between the direction of the kaon in the D° rest frame and 
the direction of the D° in the laboratory frame. In contrast, background events exhibit a 
pronounced peak at cos#* = 1. Therefore we require cos (9* < 0.9. 

In the 5-prong decay mode, the procedure is similar but, due to higher combinatorial back- 
ground, some cuts are tightened. To form the D° candidate, four tracks are combined if the 
charges of the tracks add up to zero. One track should be identified kaon, i.e. p£ B/d * > 10%, 
and this track's dE/dx probability J> dE / Ax for the pion hypothesis should be less than 10%. 
For the other three tracks, j>^ E / dx j g re q Urrec [ to be larger than 0.5%. If the mass of the D° 
candidate lies in the range 

1830 MeV < Mf^, nd < 1900 MeV, (8) 

a fifth track is added, with V dE ^ dx > 0.5% and a charge opposite to the charge of the kaon 
candidate, to form the D* candidate. 

To further reduce the combinatorial background in both decay modes, a minimum transverse 
momentum pJjT of the D* of 2 GeV is required. To ensure that the tracks forming the D* 
candidates are mostly contained in the tracking chambers, the pseudorapidity rp of the D* is 
required to be within \rp \ < 1.5, with rp = — lntan(#/2). The angle 9 is the polar angle of 
the D* candidate. 

In about 30% (8%) of the events with 5-prong (3-prong) candidates, more than one D* 
candidate passes the above cuts on and rp* . Since the probability to correctly reconstruct 
two different D* mesons in one event is negligibly small, only the D* candidates with Mp a nd 
closest to the D° mass of 1864.6 MeV p6| are retained in events with more than one D* 
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candidate. It has been checked that this method does not produce any biases. In approximately 
6% of the events, two or more D* candidates in an event share the same D° candidate, but 
different tracks were assigned as slow-pion candidate. All of these D* candidates are kept. 

Fig. |l] shows the difference between the D* and the D° candidate mass for both decay 
channels for events with AM = Mgl nd — M£jo nd < 200.5 MeV. A clear peak is observed around 
AM = 145.4 MeV which is the mass difference between the D* and the D° meson 
a background function, 

/(AM) = a • (AM - m„)\ (9) 

is performed to the upper sideband of the signal, defined by 160.5 MeV < AM < 200.5 MeV, 
where m n is the pion mass and a and b are free parameters. The \ 2 °f the fit is 13 for 18 
degrees of freedom. The fit result is superimposed for the whole AM range. In the signal 
region, defined as 142.5 MeV < AM < 148.5 MeV, a number of 100.4 ± 12.6 (stat) D* mesons 
is obtained after subtracting the fitted background from the total number of events in the signal 
region. The distribution of the wrong-charge background is also shown. It is obtained from 
the data applying identical cuts as for the signal, but requiring that the charges of the tracks 
forming the D° candidate should add up to —2 instead of 0. In addition, in the 5-prong mode 
the three pion tracks should not have equal charges. In the upper sideband, the wrong-charge 
sample gives a good description of the shape and normalisation of the background in the signal 
sample. Hence, no normalisation is applied to the wrong-charge sample. 

Using the Monte Carlo simulations, D* mesons produced in e + e~ annihilation events are 
found to contribute only around 1% to the D* signal. It was checked that the e + e~ background 
is negligible for all values of p®* ■ Non-photon-photon four-fermion background is also found to 
be negligible. 



26|. A fit of 



6 Separation of direct and single-resolved events 

We use two different methods to study the relative contributions of the direct and single-resolved 
processes to the data sample. First, we study di-jet events using the method described in more 



detail in Ref. | 2"7fl . In di-jet events, two experimental variables can be defined, x^ and x~, which 
are measures of the photon momenta participating in the hard interaction. They are calculated 
using the relation 

= w.) (10) 



where p z is the momentum component along the z axis of the detector and E is the energy of 
the jets or hadrons, respectively. Assuming in the LO picture that the two jets contain all the 
decay products of the two charm quarks, we expect for direct events that the whole energy of 
the event is contained in the two jets, i.e. x^ and x~ are close to 1. In resolved events, there 
is also energy outside the two jet cones due to the photon remnant (s). Events where either x+ 
or x~ is much smaller than 1 are expected to originate from single-resolved processes. Events 
where both x+ and x~ are much smaller than 1 are expected to originate from double-resolved 
processes. The validity of this expectation has been demonstrated in Ref. []27|] . 

The second method can be used for all D* event, not just the di-jet sub-sample. We can 
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reconstruct the scaled D* transverse momentum x^* which is given by 

*?- = ^. (id 

VIS 

If p®* is a good estimate of the transverse momentum of the charm quarks and if the charm 
quarks are produced centrally (rj = 0), the variable x®* is equal to x±. This variable is therefore 
sensitive to the ratio of the direct and the single-resolved process. As in the case of the x^ 
distribution, the direct contribution dominates at high values of x° , whereas the single-resolved 
events are concentrated at small predicted by the Monte Carlo. 

In order to reconstruct jets, a cone jet finding algorithm is applied to the signal events. As 
in the calculation of W v i S , the energy and positions of all clusters in the ECAL, the HCAL, the 
FD and the S W calorimeters and the momenta of all tracks are used in the jet finding after 
applying the matching algorithm [^4| to avoid double counting of particle momenta. The cone 



size R— J (A77) 2 + (A0) 2 is set equal to 1, where r\ and denote the pseudorapidity and the 

azimuthal angle, respectively. The minimum transverse jet energy E^ is required to be greater 
than 3 GeV. The pseudorapidity of the reconstructed jets must be less than 2. 

In Fig. 0, the fraction of events with different number of jets, n^ et , is shown for events in the 
signal region for data and for the direct and the single-resolved PYTHIA Monte Carlo samples, 
separately. About 1/3 of the data events are di-jet events. The number of signal events in 
the direct (single-resolved) Monte Carlo sample is approximately 6.5 (3.5) times the number 
of signal events in the data. In the data, the combinatorial background has been subtracted 
using the n iet distribution from the upper AM sideband. The agreement between the n iCt 
distributions in the data and the Monte Carlo is satisfactory. There are slightly more data 
events in the 3-jet bin and less data events in the 2-jet bin compared to PYTHIA. Only a small 
difference between the number of jets found in direct and single-resolved events is expected 
according to the Monte Carlo. 

Fig. ^ shows the distribution of x™ m = min(x+,x~). In the data, the combinatorial back- 
ground has been subtracted using events from the AM sideband. At small values of x™ m , the 
distribution is approximately flat, but for large x™ m values a clear enhancement is visible. This 
is expected to be due to the direct process. The PYTHIA Monte Carlo predicts 87% of the 
direct di-jet events to have x™ m > 0.7 and 82% of the single- resolved di-jet events to have 
xf n < 0.7. 

The ratio of direct to single-resolved contributions in the di-jet data is determined by a fit 
to the x™ in distribution using the method of least squares (x 2 fit). In the fit, the sum of the 
direct and single-resolved Monte Carlo samples is fixed to the number of di-jet events in the 
data, but the ratio of direct to single-resolved events is left free. According to the PYTHIA 
Monte Carlo, a significant contribution of double-resolved events in the data should show up 
as a clear enhancement at small x™ m values. As expected from Ref. ||, this is not observed, 
and the contribution from double-resolved events is therefore neglected. The fit yields that 
(46 ± 11)% of the di-jet events in the data are due to the direct and (54 ± 11)% are due to the 
single-resolved process. In Fig. || the fitted direct and single-resolved contributions are shown. 
The x 2 of the fit result is 6.1 for 9 degrees of freedom, and the fit result gives a good description 
of the data. At high x™ in values, where the direct events are concentrated, the data seem to be 
slightly shifted towards smaller x™ m values compared to the Monte Carlo. 
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In Fig. [|, the scaled D* transverse momentum x°* is plotted for all signal events after 
subtracting the combinatorial background. The ratio of direct to single-resolved contributions 
in the data is again determined by a fit using the same procedure as applied to the x™ m 
distribution. Poisson errors are assigned to the bins without data entry. The fit yields that 
(51 ±9)% of the signal events in the data are due to the direct and (49 ±9)% are due to the 
single-resolved process, consistent with the result of the fit to the a;™ m distribution. In Fig. [| 
the fitted direct and single-resolved contributions are shown. The \ 2 °f the fit result is 6.5 for 
13 degrees of freedom. Again, the fit result gives a good description of the data. 

Since in the fit to the xijT distribution all signal events are used, whereas in case of the x™ m 
distribution the fit is applied only to a part of the signal events, we use the result of the fit to 
the xljT distribution in the further analysis. 

Table [I] summarises the D* selection efficiencies e for the 3-prong and the 5-prong decay 
modes in direct and single-resolved events. The efficiencies are calculated using the Monte 
Carlo by dividing the number of reconstructed D* mesons by the number of generated D* 
mesons with p®* > 2 GeV and |?7 D | < 1.5. In the 5-prong mode, the selection efficiency for 
direct events is slightly higher than for single-resolved events, whereas in the 3-prong mode 
the selection efficiencies are about equal for direct and single-resolved event. The fractions of 
direct and single-resolved e + e~— > e + e~D*X events in the kinematical region p^ > 2 GeV and 
It? I < 1.5 are therefore assumed to be unchanged by the efficiency correction. To determine the 
systematic uncertainty, half the difference between the direct and single-resolved efficiencies in 
the 5-prong mode is used and no systematic uncertainty is used for the 3-prong mode. Together 
with the relative rate of D* mesons decaying in the 5-prong mode and the 3-prong mode, this 
yields a 8% relative uncertainty which is added quadratically to the 18% relative error of the 
direct and single-resolved contributions determined by the fit to the x^* distribution. The 
direct contribution of the process e + e~— > e + e~D*X in the kinematical region p®* > 2 GeV and 
|t7 D * I < 1.5 is determined to be rdir = (51 ± 10)% and the single-resolved process contributes to 
a fraction of 1 — r^ T = (49 ± 10)%. No significant double-resolved contribution is observed. 



7 Differential D* cross-sections 



We determine the differential cross-sections da/dp^ and dcr/d|?7 D *| for the production of D* 
mesons in anti-tagged e + e~— > e + e~D*V events as a function of the transverse momentum p®* 
and the pseudorapidity \q |. Table ^summarises the background-subtracted number of 
D* mesons and the differential cross-section da/dp^ for both decay modes. At large p% the 
statistical errors are large for the upper AM sidebands of the signal and for the wrong-charge 
distributions. These two distributions are therefore combined to determine the background in 
each decay mode and p®* bin. The background function /(AM) = a ■ (AM — m n ) b is fitted to 
this upper sideband distribution and the number of background events is calculated from the 
fit result. 

For each decay mode and for each bin in p®* , the differential cross-section da/dp^f is 
calculated using the relation 

Hn- /V rec 

d<y vv D , (12) 



dp?* e • BR • £ • Ap®* ' 
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The efficiency e is determined using the Monte Carlo by fixing the ratio of direct to single- 
resolved events to the result obtained in the previous section. The branching ratios BR(D* + — > 
K-7T+7T+) = 0.02630 ± 0.00082 and BR(D*+ -> K-7r+7r-7r+7r+) = 0.0519 ± 0.0029 are taken 
from Ref. |26j , £ is the total integrated luminosity and Ap^ is the width of the p®* bin. The 
results of both decay modes agree within the statistical uncertainties. The combined differential 
cross-section da/dp^ is given in Table The average transverse momentum (p® ) for each 
bin is determined using the method proposed in Ref. [28] . 



In Fig. [5], the combined differential cross-section da/dp^ is compared to the NLO calcula- 
tion by Frixione et al. || using the massive approach and to the NLO calculation by Kniehl 
et al. |29| using the massless approach, which was repeated by the authors specifically for the 



kinematical conditions of this analysis. In both calculations, the charm quark mass is taken 
to be m c = 1.5 GeV and the charm fragmentation is parametrised by the Peterson fragmen- 
tation function. The Peterson fragmentation parameter e c and the fraction /(c — > D* + ) of 
charm quarks fragmenting into D* + meson are e c = 0.116, f(c — > D* + ) = 0.267 in the massless 
calculation and e c = 0.035, /(c — > D* + ) = 0.233 in the massive calculation. For the massless 
calculation, the parameters were determined via a NLO fit [J2I1 to LEP1 data on D* production 
in e + e~ annihilation measured by OPAL ||30|| . The renormalisation scale /zr and the factorisa- 



tion scale jip are in both calculations defined as //r = /if/2 = £mx with mx = Jp\ + ml and 



£ = 1, where Pt is the transverse momentum of the charm quark. The GRV |2l| parametrisation 
of the parton distributions of the photon is used in the massless calculation and the GRS [[H[ 
parametrisation in the massive calculation. Despite the low transverse momenta studied, the 
agreement between data and the massless calculation is good. The massive calculation agrees 
with the data cross-section for pJjT > 3 GeV, but underestimates the data in the region of 
small p®* . The scale dependence on d<r/dp°* determined by using £ = 1/2 and £ = 2 is ap- 
proximately 10% for both calculations. The corresponding curves for the massless case are also 
shown in Fig. [5]. In addition, the massless calculation was performed using the AFG [^ and 



GS [33] parametrisations. In the massive calculation AFG and GRV were used as alternative 



parametrisations. The change of the cross-section is approximately 10% in both calculations. 

In Table ^ the number of reconstructed D* mesons with 2 GeV < pj.* < 12 GeV and 
|r] D *| < 1.5 and the corresponding differential cross-sections d<7/d|7? D | are given as a function of 
\t] D *\ for both decay modes. The numbers are determined in the same way as described above 
in the case of p^ . Within the statistical uncertainties, both decay modes yield comparable 
results. The combined differential cross-section is also given in Table |3| and plotted in Fig. |^. 
The distribution is dominated by the events at low p^* , and within the error, it is independent of 
|?7 D *|. The centres of the bins are taken as the average (|?7 D |) values. The massless calculation 
by Kniehl et al. is in good agreement with the measured differential cross-section, whereas the 
massive calculation of Frixione et al. underestimates the data, as seen already in Fig. |^. For 
the massive calculation, two additional curves are shown in Fig. ^] representing different charm 
quark masses m c with renormalisation scales /zr and factorisation scales /if as indicated in the 
figure. The combination of a small charm quark mass (m c = 1.2 GeV) with a special choice 
of the renormalisation scale (//r = 2m^ for the direct and /xr = tot/2 for the single- resolved 
process) yields a cross-section which is closer to the data, but still slightly low. 

For the determination of the systematic uncertainties, each decay mode and each bin in p®* 
or |// D *| is treated individually. The following errors are taken into account: 



12 



The relative uncertainties on BR(D* + — > K and BR(D* + — > K 7r + 7r 7r + 7r + ) of 

3.1% and 5.6%, respectively |2f 



The relative uncertainty on the selection efficiencies due to the limited number of Monte 
Carlo events varies between 5% in the lowest p^ bin for the 3-prong decay mode and 
13% in the highest pJjT bin for the 5-prong decay mode. The corresponding errors for the 
|?7 D *| distribution are 5% to 7% in all \r] D *\ bins except for the bin 1 < \i] D *\ < 1.5 in the 
5-prong mode where the error is 12%. 

The uncertainty on the number of background events determined from the fit of the back- 
ground function to the sum of the sidebands of the signal data and of the wrong-charge 
distribution. A modified background function is constructed, defined by the requirement 
Xm'm i s ^ e minimum x 2 of the fit. The relative difference between the 
number of background events determined with the modified background function and the 
number of background events determined with the original background function is taken 
as the error. Depending on the decay mode and on the p®* or \rf | bins, the relative 
uncertainties vary between 5% and 25%. 

The contributions of the direct and single-resolved Monte Carlo samples have been varied 
between 40% and 60%. In the 3-prong and in the 5-prong mode, the corresponding errors 
on the cross-sections are smaller than 7% for all bins in p®* and in |?7 D *|. 

Uncertainties in the modelling of the tracking in the central detector are assessed by 
repeating the analysis with the tracking resolutions varied in the Monte Carlo by ±10% 
around the values that describe the data best. The efficiencies obtained are compared 
with the original values, and the relative difference is quoted as the systematic error. 
Depending on decay mode and bin, this error lies between 5% and 15%. 

Uncertainties in the dE/dx probabilities for identifying kaons. In a former OPAL analysis 
|34| , D* mesons are reconstructed in the 3-prong mode using a similar set of cuts to this 



analysis. The relative error on the dE/dx probability for identifying kaons is determined 
to be around 3%. In the 5-prong mode, not studied in Ref. f34|, the pion tracks are also 



identified using the dE/dx probabilities, so the corresponding uncertainty is assumed to 
be 5%. 



The uncertainty on the integrated luminosity L is smaller than 1% and is therefore not taken 
into account. The individual systematic uncertainties are added quadratically, separately for 
each decay mode as well as for each bin in p^f or |^ D *|. The integrated cross-section cx meas of 
the process e + e~— > e + e~D*X in the kinematical region 2 GeV < p®* < 12 GeV and |r/ D *| < 1.5 
is determined to be cr°* as = 29.4 ± 3.4(stat) ± 2.4(sys) pb. 

The LO cross-sections for the direct process, cxJJ, and for the single-resolved process, cr°*, 
calculated with the PYTHIA Monte Carlo for different LO parametrisations of the parton 
densities (SaS-lD [T|, GRV JST] and LAC1 [|35j) are given in Table |. In the PYTHIA Monte 



Carlo, the charm quark mass m c was varied between 1.3 and 1.7 GeV. Since the ratio of 
direct to single-resolved cross-sections is about 1:1 in the data, the direct cross-section is well 
described by PYTHIA and the single-resolved cross-section is best described using GRV. The 
single-resolved cross-section is underestimated using SaS-lD, and the LAC1 parametrisation 
overestimates the single-resolved cross-section. 
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8 Total cross-section cr(e + e — ► e + e cc) 

For the determination of the total cross-section <r(e + e _ — > e + e~D*X), the Monte Carlo is used 
to extrapolate to the full kinematical region using the relation 

a(e + e" -> e + e"D*X) = a(e + e" -> e + e-D*X) dir + a(e + e" -> e + e"D*X) rcs 

= ^ mC as(r d i r -^ r c + (l-r dir )-^ s c ). (13) 

where 7^° and are the extrapolation factors. This allows the total cross-section of the 
process e + e~— > e + e~cc for Q\ < 4.5 GeV 2 to be calculated using the equation 

a(e + e - - e + e-cc) = 2 _ \ p ^ +) • a(e + e" - e + e"D*X). (14) 



In a previous publication J3§, the product P c = f(c -> D*+) x BR(D* + -> D°tt + ) x BR(D° -> 
K~7T + ) was derived from measurements of D* production in e + e~ collisions at v^ee = 10.5 
and 30 GeV to be P c = (7.1 ± 0.5) • 10 -3 . With the branching ratios taken from Ref. [26|, a 



hadronisation fraction /(c — > D* + ) = 0.270 ± 0.019 ± 0.010 is derived, where the last error is 
due to the branching ratio uncertainties. Since the invariant mass range of the photon-photon 
system studied in this analysis is of the same order of magnitude as the e + e~energies mentioned 
above, this value of f(c — > D*) is used in the analysis. 

The extrapolation factors for the direct events and for the single-resolved events 
are defined as the ratio of the number of all generated D* mesons in the full kinematic range 
of an d \v D *\ divided by the number of generated D* mesons with 2 GeV < < 12 GeV 
and \r] B | < 1.5. The extrapolation factors are TZ^f = 12.6 and TZ%& = 18.4 obtained using 
the combination m c = 1.5 GeV and the Peterson fragmentation with e c = 0.031. 

The extrapolation introduces systematic uncertainties due to the modelling of the fragmen- 
tation of the charm quarks into D* mesons which influence mainly the distributions. To 
determine the systematic errors on the extrapolation factors, different m c values and different 
fragmentation functions were used for the event generation in the Monte Carlo: 

• The charm quark mass m c was varied between 1.3 and 1.7 GeV. 

• e c = 0.0851 |29| was used in the Peterson fragmentation function. 



The Lund symmetric fragmentation function [37] was used with the parameters a = 1.95 



and b = 1.58 determined in Ref. [35]. 



In case of the single-resolved process, the GRV parametrisation was used as an alternative 
parametrisation in combination with different charm quark masses and fragmentation 
functions. 



For all studied combinations of m c , fragmentation functions and parametrisations of the 
parton densities, the direct and the single-resolved Monte Carlo samples were added in such a 
way that in the kinematical range 2 GeV < < 12 GeV and \rp \ < 1.5 the cross-section 
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is equal to cx meas , the direct contribution is r dir = 51% and the single-resolved contribution is 
I - r dir = 49%. The resulting differential cross-sections da/dp^* and dcr/d|?7 D *| are in agreement 
with the measured differential cross-sections (Figs. [5|-§|). Therefore the mean quadratic devia- 
tion of all calculated extrapolation factors from the central values 7£^ r c = 12.6 and 7^ s c = 18.4 
is used to determine the relative systematic uncertainties of 19% for 1Z$~ and 27% for TZf 10 



^rcs 



Combining Eqs. flHf ) and ([14]) , we determine the total cross-section of the process e + e — > 
e+e"cc to be a(e+e" -> e + e"cc) = 842 ± 97 (stat) ± 75 (sys) ± 196 (extr) pb at ^ ee = 
183 — 189 GeV. The first error is the statistical, the second error is the systematic error and 
the third error is the extrapolation uncertainty. The direct contribution is determined to be 
a(e + e~ — > e + e~cc)dir = 351 ± 40 (stat) ± 79 (sys) ± 66 (extr) pb and the single-resolved 
contribution to be cr(e + e~ — > e + e~cc) res = 491 ± 56 (stat) ± 111 (sys) ± 130 (extr) pb. 

The separation of direct and resolved events in heavy quark production is scheme-dependent 
in the NLO massless calculation, but it is unambiguous in LO and in the NLO massive calcula- 
tion. Using m c = 1.5 GeV, the LO direct cross-section in PYTHIA lies in the range 

300i£ pb. 

The LO calculation of Ref. @ gives 382I94 6 pb, and the NLO calculation 593tlH pb- The upper 
and lower error correspond to m c = 1.3 and 1.7 GeV, respectively. The OPAL anti-tagging 
condition was applied to these calculations. The measured direct cross-section agrees well with 
the LO calculations, whereas it lies at the lower end of the NLO calculation. This could be due 
to the separation procedure for direct and single-resolved events which uses distributions from 
a LO Monte Carlo. 

Fig. |7| shows the total cross-section <r(e + e~ — > e + e~cc) compared to other measurements 
and to the NLO calculation of Ref. [@ using the GRS parton distributions. The calculation is in 
good agreement with the OPAL result within the large band of uncertainties due to variations 
of m c , /ir, and /if- The anti-tagging condition used in this paper has been applied to the NLO 
calculation. It should be noted that the anti-tagging conditions of the different experiments 
are not identical. The OPAL result is consistent with the L3 result at v^ce = 167 GeV and 
about 1.5 standard deviations below the L3 result at \/s ee = 183 GeV fl39 |. 



9 Determination of F^ c 

In this section, deep inelastic electron-photon scattering is studied using single-tagged events. 
The D* production cross-section, as well as the charm production cross-section and the charm 
structure function K7 C of the photon are determined from events with a beam electron scattered 
into the forward detectors (tagged events). An event is tagged (SW-tagged or FD-tagged) if 
the energy of the scattered electron E', measured in the angular range 33 mrad < 6' < 55 mrad 
for the SW or 60 mrad < 9' < 120 mrad for the FD, exceeds 50 GeV in one hemisphere of 
the detector. The corresponding approximate ranges in Q 2 are 5 GeV 2 < Q 2 < 30 GeV 2 and 
30 GeV 2 < Q 2 < 100 GeV 2 for SW-tagged and FD-tagged events, respectively. The selection 
of D* candidates is identical to the selection in the anti-tagged case, with three exceptions: 



• For the calculation of the visible invariant mass W v \ s , clusters in the SW or FD in the 
hemisphere of the tag are excluded. 
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• The combinatorial background in tagged events is smaller due to the slightly smaller 
mean number of tracks per event. This makes it possible to include D* mesons with 

> 1 GeV for SW-tagged events. 

• In FD-tagged events, the D* mesons have higher transverse momenta p®* due to the 
transverse momentum balance between the tagged electron and the hadronic system. To 
improve the signal to background ratio, a cut p^* > 3 GeV is applied for FD-tagged 
events. 

Fig. [I] shows the distribution of the difference between the D* and the D° candidate mass found 
in the tagged sample. The fit of the background function Eq. || to the upper sideband of the 
signal was performed in the range 154.5 MeV < AM < 200.5 MeV. The x 2 °f the fit result is 
25 for 21 degrees of freedom. Subtracting the background predicted by the fit, 29.8 ± 5.9 (stat) 
D* mesons are found in the signal region of the tagged events. The combinatorial background 
in the upper sideband is also well described by the tagged wrong-charge sample. Background 
subtraction with the wrong-charge sample gives a consistent result for the number of D* events. 
Due to the small number of D* mesons both D* decay modes are combined for the further 
analysis. No double-tagged D* event has been found, i.e. an event with energy deposits of 
more than 50 GeV in the forward calorimeters in both hemispheres. 

Fig. H shows the distributions of W v - ls and of the measured Q 2 of the tagged signal events, 
and the p^* distribution is presented in Fig. [| The data are compared to the predictions of the 
HERWIG and Vermaseren Monte Carlo generators, normalised to the number of data events. 
Both Monte Carlo generators give a good description of the shape of the data distributions. 

We determine the cross-section for D* production in deep inelastic electron-photon scattering 
in the well-measured kinematic range: pifT > 1 GeV for an electron scattering angle 33 mrad < 
6' < 55 mrad (SW) or > 3 GeV for 60 mrad < 6' < 120 mrad (FD), \r] D *\ < 1.5 and 
E' > 50 GeV, 5 GeV 2 < Q 2 < 100 GeV 2 ; using almost the whole accessible Q range defined 
by 6' and E' . 

The analysis is performed in two bins of x with 0.0014 < x < 0.1 and 0.1 < x < 0.87. 
The x range is limited by the Q 2 range, by the minimum kinematically allowed invariant mass 
W > 3.88 GeV needed to produce a D* meson, and by the event selection cut W v i S < 60 GeV. 
To take into account the detector acceptance and resolution in x the data are corrected using a 
2x2 matrix. The measured x vis is calculated from Eq. |] using W vis and the measured value of 
Q 2 . The resolution effects in the Q 2 reconstruction are small compared to the resolution effects 
in measuring x. They can therefore be neglected. 

The D* selection efficiency for x > 0.1 is given by the ratio of the number of reconstructed D* 
mesons originating from events with x > 0.1 to all generated D* mesons in events with x > 0.1, 
in the restricted kinematic range defined above. The selection efficiency for x > 0.1 is about 
(21 ± 2)% (not including the branching ratios). For x > 0.1 the selection efficiencies obtained 
from both Monte Carlo generators are consistent, but for x < 0.1 the selection efficiencies are 
around (30 ±3)% according to HERWIG and around (18 ±2)% according to Vermaseren. Both 
programs predict that about one third of the selected D* events generated with x < 0.1 are 
reconstructed with x v - ls > 0.1, whereas migration from x > 0.1 to x v - ls < 0.1 is very small. 



16 



Table |5| summarises the number of reconstructed D* mesons and gives the measured values 
of the cross-section 

which is the deep inelastic electron-photon scattering cross-section for D* production in the 
restricted kinematic range as defined above. It is calculated from the number of D* events, 
N§?, obtained from the 2x2 matrix correction using both HERWIG and Vermaseren. For 
the combined D* branching ratios into the 3-prong and 5-prong mode, we use BR = 0.0782 ± 



0.0030 [pql . The total integrated luminosity C is 220 pb . The average of the cross-sections 



corrected with HERWIG and Vermaseren is also given in Table |5]. 

For x > 0.1, both Monte Carlo models yield consistent results. For x < 0.1, the difference 
between the cross-sections <r°* obtained using HERWIG and Vermaseren is due to the different 
D* selection efficiencies. The following systematic errors on cr^* are taken into account: 



The limited number of Monte Carlo events leads to an uncertainty of approximately 15% 

D* 
tag 



on a^* for each of the Monte Carlo generators. 



Within the statistical uncertainties, the HERWIG and Vermaseren models yield consistent 
corrected numbers of D* mesons for x > 0.1, whereas for x < 0.1, the corrected numbers 
of D* mesons obtained with the HERWIG and Vermaseren Monte Carlo models differ by 
more than 2 standard deviations. Therefore only for x < 0.1, half the difference between 
<t°* using HERWIG and Vermaseren is taken as error on the averaged value of er°* . 

The combined relative uncertainty on the branching ratios BR(D* + — > K~7r + 7r + ) and 



BR(D*+ K-7r + 7r-7r+7r+) is 3.8% |26 



The uncertainty in the number of background events in the signal region estimated in the 
same way as for anti-tagged events gives relative errors of 12% for x < 0.1 and 6% for 
x > 0.1. 

The relative uncertainty due to the modelling of the tracking in the central detector is 
estimated to be 8% using the corresponding errors determined for the differential cross- 
sections dcr/dpJjT an d dcr/d|?7 D *| for the anti-tagged events. 

The relative uncertainty due to the use of the dE/dx probabilities for identifying kaons 
and pions is estimated to be 4%, also using the corresponding errors determined for the 
differential cross-sections da/dp^* and dcr / d|?7 D * | for the anti-tagged events. 

The uncertainty due to the measurement of the energy E' of the tagged electron is assessed 
by shifting the reconstructed quantity in the Monte Carlo according to its resolution and 
by repeating the analysis. The change on <r°* is around 3% and is taken into account 
as relative error. The uncertainty due to the measurement of the scattering angle 9' is 
determined in the same way as for E' . The relative change on o"° g is found to be only 
around 1%. This error is therefore neglected. The uncertainty due to the measurement 
of the visible invariant mass W v - ls of the event is estimated to be only around 1% and is 
therefore also neglected. 
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All systematic errors are added in quadrature. 



For the determination of the total cross-section of D* production in deep inelastic electron- 
photon scattering, <r(e + e~ — > e + e~D*X), the Monte Carlo models are used to extrapolate to 
the whole kinematic region. This allows the total charm cross-section in deep inelastic electron- 
photon scattering to be calculated via the relation 



<r e e — > e e cc 



2-/(c^D*+; 

1 

2-/(c^D*+; 



• a(e + e" -> e + e"D*X) 

• « " <C (16) 



The extrapolation factor TZ^g is defined in the same way as in the anti-tagged case. Table |6| gives 
the values of the total charm cross-section, <r(e + e~ — > e + e~cc), extrapolated using HERWIG and 
Vermaseren as well as the averaged cross-section. The extrapolation error has been determined 
in the following way: 

For x > 0.1, both Monte Carlo generators predict very similar values for TZ^S (4.8/4.6 for 
HERWIG/ Vermaseren) and thus for <r(e + e - — > e + e~cc). The uncertainty on is determined 
in the same way as in the anti-tagged case. It is found that the influence of the charm quark 
mass and fragmentation function is small, and the relative uncertainty on T^Jfg and thus on 
the averaged cross-section a(e + e~ — > e + e~cc) is only 5%. 

In contrast, for x < 0.1, the Monte Carlo generators predict very different extrapolation 
factors due to the large discrepancy between the predicted invisible part of the cross-section. 
For x < 0.1, the HERWIG extrapolation factor, TZ^ = 12.9, is more than twice as large as the 
Vermaseren factor TZf^ = 5.1. The predicted cross-sections of the Monte Carlo models and the 
NLO calculation of Laenen et al. @ are given in Table 0. Since the hadron-like contribution 
is neglected in the QPM, the Vermaseren cross-section is much smaller than the LO and the 
NLO cross-section for x < 0.1. In contrast, mainly due to the massless approach taken, the 
prediction from HERWIG is higher than the cross-section from the LO and the NLO calculation. 
Therefore it is likely that the correct cross- sect ion, and therefore the correct extrapolation 
factor, lies within the range of the two Monte Carlo predictions. Half the difference between 
the two extrapolated cross-sections is taken into account as extrapolation error on the averaged 
cross-section <r(e + e~ — > e + e~cc). 

Finally, the value of the charm structure function F% c {x, (Q 2 )) of the photon, averaged over 
the corresponding bin in x, is determined by 

FUx, <Q 2 » = a(e+e- - e+e'cc) • f ^ ^ ) , (17) 

\a{e+e -> e+e cc)/ NLO 

where the ratio (F^ c (x, (Q 2 ))/a(e + e~ — > e + e~cc))Nix> is given by the NLO calculation of Laenen 
et al. ||. The mean virtuality in the measured region 5 GeV 2 < Q 2 < 100 GeV 2 is about 
(Q 2 ) ~ 20 GeV 2 , in agreement with the values from the generated HERWIG and Vermaseren 
Monte Carlo events. The F^ c {x, (Q 2 )) values are given in Table ^. They are calculated from 
the individual charm cross-sections obtained using the HERWIG and Vermaseren models and 
from the averaged cross-section. 
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In Fig. |10|a), the measured cross-sections obtained using the individual Monte Carlo models, 
are compared to the calculation of Laenen et al. |3| performed in LO and NLO and to the Monte 
Carlo results, and Fig. [10] b) shows the charm structure function F^ c . The NLO prediction 
is based on m c = 1.5 GeV and the renormalisation and factorisation scales are chosen to be 

= = Q- The calculation is obtained for the sum of the point-like and hadron-like 
contributions to F% c , using the GRV-NLO parametrisation in the calculation of the hadron-like 
part. The NLO corrections are predicted to be small for the whole x range. The NLO calculation 
is shown as a band representing the uncertainty of the theoretical prediction, evaluated by 
varying the charm quark mass between 1.3 and 1.7 GeV and by changing the renormalisation 
and factorisation scales in the range Q/2 < /ir = fip < 2 Q. 

For x > 0.1, all cross-section predictions in Fig. |10| a) are consistent with one another. 
Because the cross-section prediction from the Vermaseren model is consistent with the point- 
like contribution to the LO calculation for the whole x range, the contributions from longitudinal 
photons are expected to be small. For x < 0.1, the situation is different. The NLO calculation 
predicts the hadron-like and point-like component to be of about equal size. Therefore the 
purely point-like QPM prediction of the Vermaseren model is expected to underestimate the 
data if a hadron-like contribution exists. The HERWIG Monte Carlo predicts the highest cross- 
section, which is expected, since the massless approach should overestimate the cross-section, 
as explained in Section 4. 

The different behaviour of the Monte Carlo cross- sect ions in the two regions of x is reflected 
in the measured cross-sections shown in Fig. [10] a). For x > 0.1 the individual measured cross- 
sections obtained by correcting with HERWIG and Vermaseren are very similar, the error of 
the measured cross-section is dominated by the statistical uncertainty, and the NLO calculation 
is in good agreement with the data. In contrast, for x < 0.1, the result suffers from the strong 
model dependence discussed above. The result based on the HERWIG generator is much higher 
than the result obtained using the Vermaseren model. Despite this uncertainty the corrected 
data suggest a cross-section which is above the purely point-like component, i.e. the hadron- 
like component of Fr[ c is non-zero. This observation is independent of the Monte Carlo model 
chosen for correction. Averaging the individual results is therefore safe for x > 0.1, but for 
x < 0.1 the averaged result suffers from large model uncertainties and has to be interpreted 
with care. In Fig. [H^ b) the cross-section measurements are converted into the measured charm 
structure function using Eq. [I7|. The conclusions derived from F^ c and from the cross-sections 
are the same. 



In Fig. [11], the averaged results are presented. Fig. |] 1] a) shows the cross-section on a 



linear scale in x in comparison to the same predictions as in Fig. |10[ In Fig. |ll| b) the charm 



structure function is presented on a logarithmic scale in x for (Q 2 ) = 20 GeV 2 . The data points 
for F^ c are located at the mean value of x, denoted with (x). The values are the averaged (x) 
values obtained with both Monte Carlo generators, and half the difference of the HERWIG 
and Vermaseren predictions is taken as the uncertainty. For x > 0.1, the predicted (x) is 
around 0.32 and the difference between the HERWIG and Vermaseren programs is invisible. 
The point-like contribution decreases for decreasing x, whereas the hadron-like component 
rises. Consequently, for x < 0.1, the HERWIG Monte Carlo predicts a smaller average value of 
(x) = 0.028 than the Vermaseren Monte Carlo which yields (x) = 0.054. 
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In addition to the full NLO prediction, the predicted hadron-like component of F^ c is also 
shown in Fig. [IT] b). This contribution is very small for x > 0.1 and therefore in this range the 
NLO calculation is an almost purely perturbative prediction with the charm quark mass and 
the strong coupling constant as the only free parameters. This prediction nicely describes the 
data. To illustrate the shape of K7 C the data are also compared to the GRS-LO prediction 
and to the point-like component alone both shown for Q 2 = 20 GeV 2 . The point-like component 
strongly decreases for decreasing x. The full F^ c evaluated at Q 2 = 20 GeV 2 agrees with the 
data. The change of within the range of Q 2 studied is large. The maximum value of F^ c 
for x > 0.1 rises by about a factor of five between Q 2 = 5 GeV 2 and Q 2 = 100 GeV 2 and the 
charm thresholds moves from about x = 0.35 to about x = 0.9. 

In conclusion, for x > 0.1, the purely perturbative NLO calculation is in good agreement 
with the measurement and for x < 0.1, the measurement suffers from large uncertainties of the 
invisible cross-section predicted by the HERWIG and Vermaseren Monte Carlo models, and 
therefore the result is not very precise. However, despite the large error in this region the data 
suggest a non-zero hadron-like component of K7 C . 



10 Conclusion 



We have measured the inclusive production of D* 1 * 1 mesons in photon-photon collisions using 
the OPAL detector at LEP at e + e~ centre-of-mass energies \/s ee = 183 and 189 GeV. The D* + 
mesons are reconstructed in their decay to D°7r + with the D° observed in the two decay modes 
K~7r + and K~7i + n^n + . In total, 100. 4± 12.6 (stat) D* mesons are selected in anti-tagged events 
and 29.8 ± 5.9 (stat) D* mesons in single-tagged events. 

In the anti-tagged event sample, the direct and single-resolved contributions are separated 
using di-jet events reconstructed with a cone jet finding algorithm, and for all observed events 
by fitting the distribution of the scaled D* transverse momentum x^ . Both methods yield 
consistent results, and due to the larger statistics used, the second method is more precise. It is 
found that in the kinematical region > 2 GeV and |^ D * | < 1.5 the direct contribution to the 
process e + e~— > e + e~D*X is (51 ± 10)% and that the single-resolved contribution is (49 ± 10)%. 

Differential cross-sections as functions of the D* transverse momentum and pseudorapidity 
are measured for anti-tagged events and are compared to a NLO calculation by Kniehl et al. 
using the massless approach, and by Frixione et al. || using the massive approach. It is found 
that despite the low values of studied the massless calculation is in good agreement with 
the data. The massive calculation agrees with the measured cross-section for p!p* > 3 GeV but 
underestimates the data for lower values of p^* . 

The total cross-section of the process e + e~— > e + e~cc, where the charm quarks are produced 
in the collision of two quasi- real photons, is measured to be <r(e + e~ — > e + e~cc) = 842 ± 
97 (stat) ± 75 (sys) ± 196 (extr) pb, with a direct contribution of a(e + e _ — > e + e~cc)d ir = 
351 ± 40 (stat) ± 79 (sys) ± 66 (extr) pb and a single-resolved contribution of cr(e + e - — > 
e+e-cc) res = 491 ± 56 (stat) ± 111 (sys) ± 130 (extr) pb. The NLO calculation of Ref. @ and 
the measurements by L3 |39| are in agreement with this result. 
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The first measurement of the charm structure function -F^ c (a;, (Q 2 )) of the photon has been 
performed based on 29.8 ± 5.9 (stat) D* mesons reconstructed in single-tagged events. The 
value of F2 yC (x, (Q 2 )) is determined for an average (Q 2 ) of 20 GeV 2 and in two regions of x, 
0.0014 < x < 0.1 and 0.1 < x < 0.87. The NLO corrections to F^ c (x, (Q 2 )) are predicted to be 
small for all x and the contribution of the hadron-like component is negligible for x > 0.1, which 
means that F2 C (x, (Q 2 )) can be predicted purely perturbatively in this region. For x > 0.1, the 
perturbative NLO calculation of Laenen et al. [|Tj is in good agreement with the measurement. 
For x < 0.1, the measurement suffers from large uncertainties of the invisible cross-section 
predicted by the HERWIG and Vermaseren Monte Carlo models, and therefore the result is 
not very precise. However, despite the large error in this region the data suggest a non-zero 
hadron-like component of F^ c . 
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D° -> K-it+ 


D° -> K-7T+7r-7r+ 


direct 


(40.8 ± 1.7)% 


(14.4 ±0.8)% 


single-resolved 


(38.1 ±2.0)% 


(11.2 ±0.9)% 



Table 1: D* selection efficiencies e for the two decay modes and for direct and single- resolved 
events as determined from the PYTHIA Monte Carlo (for anti-tagged events only). The selec- 
tion efficiencies refer to D* mesons with pJjT > 2 GeV and \i] D *\ < 1.5. Only statistical errors 
are given. 





(P?*> 


Arrec 

jv d * 


da/dp®* [pb/GeV] 


[GeV] 


[GeV] 


K-7T + 




K-7T + 




combined 


2-3 


2.46 


42.9 ±7.7 


27 .4 ±7.3 


20.6 ±3.7 ±3.0 


21.5 ±5.8 ± 6.1 


20.8 ± 3.1 ± 2.4 


3-5 


3.82 


18.4 ±4.5 


11.4 ±4.2 


3.7 ±0.9 ±0.5 


3.1 ±1.2 ± 0.7 


3.5 ± 0.7 ± 0.3 


5-12 


7.30 


8.3 ±3.0 


4.5 ±2.2 


0.38 ±0.14 ±0.06 


0.25 ±0.12 ±0.05 


0.31 ±0.09 ±0.09 



Table 2: Number of reconstructed D* mesons with \rj \ < 1.5 in bins of p3f for both decay 
modes after background subtraction (for anti-tagged events only). The differential D* cross- 
section as a function of pJjT for each decay mode and the combined cross-section is also given. 
The first error is statistical and the second error is systematic. 



24 



\ri D *\ 


iv D* 


da/d^ *! [ P b] 


K _ 7T + 


K 7T + 7T 7T + 


K _ 7T + 


K 7T + 7T 7T + 


combined 


0.0-0.5 


29.1 ±5.9 


17.0 ±5.3 


22.6 ±4.6 ±2.7 


18.0 ± 5.6 ±4.7 


21.0 ±3.5 ±2.1 


0.5-1.0 


18.1 ±5.3 


24.0 ±5.8 


14.5 ±4.2 ±3.8 


23.8 ± 5.8 ±4.9 


18.0 ±3.4 ±2.8 


1.0- 1.5 


22.8 ±5.0 


6.3 ±3.7 


25.3 ±5.6 ±2.9 


19.5 ± 11.7 ± 6.1 


24.2 ±5.0 ±2.2 



Table 3: Number of reconstructed D* mesons with 2 GeV < p® < 12 GeV in bins of \rj | for 
both decay modes after background subtraction (for anti-tagged events only). The differential 
D* cross-section as a function of \i] D *\ for each decay mode and the combined cross-section is 
also given. The first error is statistical and the second error is systematic. 





°Z [P b ] 


°?cs [Pb] 


PYTHIA 


14.0 - 14.9 


3.8 - 5.5 
(SaS-lD) 


6.9 - 9.9 
(GRV) 


23.0 - 37.0 
(LAC1) 



Table 4: Predicted integrated LO cross-section of the process e + e — > e + e D*X in the kinemat- 
ical region 2 GeV < p?* < 12 GeV and \r] D *\ < 1.5 calculated with the PYTHIA Monte Carlo 
using different parametrisations of the parton densities (anti-tagged events only). Direct and 
single-resolved cross-sections are given separately. The charm quark mass was varied between 
1.3 and 1.7 GeV. 



X 




a tag [pb]> corrected with 


CT te g [Pb] 


HERWIG 


Vermaseren 


average 


0.0014-0.1 


9.9 ±3.6 


3.1 ± 1.1 ±0.7 


4.7 ±1.8 ±0.9 


3.9 ± 1.4 ± 1.0 


0.1-0.87 


20.0 ±4.7 


4.0 ± 1.4 ±0.7 


4.0 ±1.4 ±0.7 


4.0 ± 1.4 ±0.6 



Table 5: Number of reconstructed D* mesons, N^*, found in tagged events and a® , obtained 
by correcting with HERWIG and Vermaseren in two bins of x, where <r°* is the deep inelastic 
electron-photon scattering cross-section for D* production in the restricted kinematic range as 
defined in the text. The first error on o"° g is the statistical error of the data and the second 
error is the systematic error. The right hand column gives the averaged cross-sections obtained 
by correcting with the two Monte Carlo models. 
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X 


<r(e + e — > e + e cc) [pb], corrected with 


cr(e + e — > e + e cc) [pb] 


HERWIG 


Vermaseren 


average 


0.0014-0.1 


92.6 ±34.1 ±24.6 


36.9 ± 13.6 ±9.8 


64.8 ± 23.9 ± 17.2 ±27.9 


0.1-0.87 


35.5 ± 12.4 ± 5.2 


34.3 ± 11.9 ±5.1 


34.9 ± 12.1 ± 5.1 ± 1.7 


X 


Fq,,c{ x i {Q 2 ))/ a i corrected with 


Fl(x, (Q 2 ))/a 


HERWIG 


Vermaseren 


average 


0.0014-0.1 


0.39 ±0.14 ±0.10 


0.16 ±0.06 ±0.04 


0.27±0.10±0.07±0.12 


0.1-0.87 


0.11 ±0.04 ±0.02 


0.11 ±0.04 ±0.02 


0.11 ±0.04 ±0.02 ±0.01 



Table 6: Total charm cross-section in deep inelastic electron-photon scattering, a(e + e~ — > 
e + e~cc), for 5 GeV 2 < Q 2 < 100 GeV 2 and the charm structure function of the photon divided 
by the fine structure constant, F^ c (x, (Q 2 ))/a, averaged over the corresponding bin in x for 
(Q 2 ) = 20 GeV 2 . The cross-section and F^ c (x, (Q 2 ))/a are presented corrected using both the 
HERWIG and Vermaseren Monte Carlo models. The averaged values for the cross-section and 
K7 c (x, (Q 2 ))/a are also given. The first errors are statistical, the second errors systematic, and 
the third errors are the extrapolation uncertainties. 



X 


cr(e + e — > e + e cc) [pb] 


F 2 ] c (x, (Q 2 ))/a 


HERWIG 


Vermaseren 


LO 


NLO 


LO 


NLO 


0.0014-0.1 


22.6 


7.7 


15.3 


16.3±i;f 


0.070 


0.069j£g£ 


0.1-0.87 


20.3 


24.7 


26.1 


30.1±I;1 


0.082 


0.0971^ 



Table 7: Predicted total charm cross-section in deep inelastic electron-photon scattering, 
a(e + e~ — > e + e~cc), for 5 GeV 2 < Q 2 < 100 GeV 2 according to the Monte Carlo genera- 
tors HERWIG and Vermaseren and according to the calculation of Laenen et al. performed 
in LO and NLO. For the calculation, also the F^ c (x, {Q 2 ))/a prediction is quoted in LO and 
NLO for Q 2 = 20 GeV 2 . The errors of the NLO results are obtained by varying the charm 
quark mass and the renormalisation and factorisation scales in the calculation. 
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Figure 1: Mass difference AM = M^ nd — Mpo nd for both decay modes for the anti-tagged and 
tagged sample. In both samples, a clear peak is visible around AM = Md* — Mqo = 145.4 MeV. 
The result of a fit of the background function /(AM) = a - (AM — m 7T ) b to the upper sidebands 
is superimposed. The fit regions are AM > 160.5 MeV for the anti-tagged events and AM > 
154.5 MeV for the tagged events. The open histograms represent the corresponding wrong- 
charge background samples which give a good description of the combinatorial background. 
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Figure 2: Fraction of signal events with different n iCt determined with the cone jet finding 
algorithm (for anti-tagged events). The dots represent the data after subtraction of the com- 
binatorial background. The solid line shows the PYTHIA prediction for the direct and the 
dashed line for the single-resolved sample, respectively. 



42 18 

> 14 
12 
10 
8 
6 
4 
2 




t — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

OPAL 

-f- data 

m PYTHIA: 46% direct 
□ PYTHIA: 54% single-resolved 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



X mm 



Figure 3: Minimum of x+ or x~ for di-jet events in the signal region (anti-tagged events only). 
The data, represented by the dots, are background subtracted using sideband events. The 
enhancement at large values of a:™ 111 is due to the direct process. The histograms are the result 
of a fit of the relative contributions of the direct and single-resolved Monte Carlo samples to 
the data. The open histogram shows the single-resolved, the hatched histogram the direct 
contribution to the fit result. 
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Figure 4: Scaled transverse momentum of the D* meson, x° , for all signal events (anti- 
tagged events only). The dots represent the background subtracted data. The histograms are 
the result of a fit of the relative contributions of the direct and single-resolved Monte Carlo 
samples to the data. The open histogram shows the single-resolved, the hatched histogram the 
direct contribution to the fit result. 
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Figure 5: The differential D* cross-section, da/dp® , for the process e + e~— > e + e~D*X in the 
range |r] D *| < 1.5 (for anti-tagged events). The dots represent the combined cross-sections from 
both decay modes. The inner error bars give the statistical error and the outer error bars the 
statistical and the systematic error added in quadrature. The data are compared to a NLO 
calculation by Kniehl et al. using the massless approach for three different renormalisation 
and factorisation scales, /ir. and fip, and to a NLO calculation by Frixione et al. using the 

massive approach. The quantity wt is defined as mx = \Jvr + m 1 where pt is the transverse 
momentum of the charm quark. 
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Figure 6: The differential D* cross-section dcx/d|?7 D * | for the process e + e~— > e + e~D*X in the 
range 2 GeV < plf < 12 GeV (for anti-tagged events). The dots represent the combined cross- 
sections from both investigated decay modes. The inner error bars give the statistical error and 
the outer error bars the statistical and the systematic error added in quadrature. NLO QCD 
calculations by Kniehl et al. using the massless approach are also shown as well as NLO QCD 
calculations by Frixione et al. using the massive approach using different renormalisation scales 
separately for the direct (dir) and the single-resolved (res) contributions and different charm 

quark masses. The quantity m T is defined as m T = + m%, where is the transverse 
momentum of the charm quark. 
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Figure 7: Comparison of measured cross-sections for the process e + e~— > e + e~cc where the 
charm quarks are produced in the collision of two quasi-real photons. The outer error bars on 
the OPAL points represent the total errors, including the extrapolation uncertainty, and the 
inner bars are the statistical errors. The values for TASSO, TPC/27, JADE, TOPAZ, AMY 
and VENUS are taken from Ref. |£jg, for ALEPH from Ref. |TD[ and for L3 from Ref. 
The band shows a NLO calculation of the process e + e~— > e + e~cc [0] for a charm quark mass 
between 1.3 and 1.7 GeV using the GRS parametrisation for the parton distributions of the 
photon. 
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Figure 8: The distributions of a) the visible invariant mass, W V - 1S , and b) the negative four- 
momentum squared, Q 2 , for the tagged signal events. The data are compared to the predictions 
of the HERWIG and Vermaseren generators. The Monte Carlo distributions are normalised to 
the number of data events. 



03 

£ 12 
10 

8 



6 E- 
4 
2 




n i i i i i i i i i i i i i i i i r 



OPAL 



data 

HERWIG 
Vermaseren 



10 



12 



P [GeV] 

T 

Figure 9: The distribution of transverse momentum pJjT f° r the D* mesons in the tagged signal 
events. The data are compared to the predictions of the HERWIG and Vermaseren generators. 
The Monte Carlo distributions are normalised to the number of data events. 
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Figure 10: Results compared with predictions for a) the deep inelastic electron-photon scat- 
tering cross-section cr(e + e~ — > e + e~cc), with 5 GeV 2 < Q 2 < 100 GeV 2 and b) for the charm 
structure function of the photon divided by the fine structure constant, F$ c (x, (Q 2 ))/a at 
(Q 2 ) = 20 GeV 2 . The data are shown individually corrected with the HERWIG and Ver- 
maseren Monte Carlo generators. The inner error bar is the statistical and the outer error bar 
is the full error. The measurements are presented at the central x values of the bins. The 
results obtained with the HERWIG and Vermaseren generators are slightly separated for a 
better visibility. The calculation of Laenen et al. [[| is performed in LO and NLO. The band 
for the NLO calculation indicates the theoretical uncertainties assessed by varying the charm 
quark mass and renormalisation and factorisation scales. 



34 



^100 

> CJ 

+ a, 75 







; 


1 


<Q 2 > = 20 GeV 2 


l 

-T" 


OPAL (a) _E 










NLO (Laenen et al.) j 











LO J 


: < 


> 






HERWIG = 








■f 


Vermaseren E 








J 




- i i 


i i 1 


i i i i I i i i i I i i i i I i i i 


i I i i i i I 


i i i i l i i i i l i i i i l i i i i - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 



A 

O 0.4 

V 



0.2 - 







- 1 

: <Q 2 > = 20 GeV 2 
=_ — GRS-LO 
• • • point-like 

: l — < 


-f OPAL (b) : 

NLO (Laenen et al.) : 

- - NLO hadron-like : 

— LO = 
>— 1 : 









10 



-3 



10 



-2 



10 



-1 



Figure 11: OPAL results for a) the deep inelastic electron- photon scattering cross-section 
a(e + e~ — > e + e _ cc), with 5 GeV 2 < Q 2 < 100 GeV 2 and b) for the charm structure func- 
tion of the photon divided by the fine structure constant, F^ c (x, {Q 2 ))/a, for an average (Q 2 ) 
of 20 GeV 2 . The data points are obtained averaging the results obtained with the HERWIG 
and Vermaseren Monte Carlo models. The outer error bar is the total error and the inner 
error bar the statistical error. The x values of the data points are obtained by averaging the 
mean x values taken from the HERWIG and Vermaseren generators. The data are compared 
to the calculation of Laenen et al. M performed in LO and NLO. The band for the NLO 
calculation indicates the theoretical uncertainties assessed by varying the charm quark mass 
and renormalisation and factorisation scales. In a) the cross-section predictions of the Monte 
Carlo generators HERWIG and Vermaseren are also given, b) also shows the prediction of the 
GRS-LO parametrisation for the whole structure function at (Q 2 ) = 20 GeV 2 and the point-like 
component alone. 
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